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We present several possibilities offered by the reaction dynamics of dissipative heavy 
ion collisions to study in detail the symmetry term of the nuclear equation of state, EOS. 
In particular we discuss isospin effects on the nuclear liquid-gas phase transition, Isospin 
Distillation, and on collective flows. We stress the importance of a microscopic relativistic 
structure of the effective interaction in the isovector channel. The possibility of an early 
transition to deconfined matter in high isospin density regions is also suggested. We finally 
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select Eleven observables, in different beam energy regions, that appear rather sensitive 
to the isovector part of the nuclear EOS, in particular in more exclusive experiments. 



Introduction 

There are quite stimulating predictions on new phases of Asymmetric Nuclear Mat- 
ter, ANM, that eventually could be reached during heavy ion reaction dynamics with 
radioactive beams [ HJ 13 El IH El • More symmetric and narrower isotopic distribu- 
tion of primary fragments are predicted, and sensitive to the symmetry term of the EOS. 
For semi-central collisions the dynamics of the participant zone appears also to be quite 
affected by the symmetry term [0 El E3 ITU1 ITT] . 

Collective flows are particularly interesting since we can probe different density regions 
of the EOS. A very stimulating result shown here is the sensitivity to the microscopic 
covariant structure of the isovector channel in the in medium interaction. A related earlier 
possible transition to a mixed phase with deconfined matter is finally presented. 
The EOS symmetry term 

The behaviour of the symmetry term of the nuclear EOS is poorly known in regions 
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Figure 1. Symmetry contribution to the 
mean field at I = 0.2 for neutrons and 
protons: dashed lines "asy-soft", solid 
lines " asy-stiff" , long dashed lines "asy- 
superstiff" 



Figure 2. Density dependence of pro- 
ton and neutron chemical potential for 
asy- super stiff (solid lines) and asy-soft 
(dashed line) EOS. 



far from normal density. In the following we will compare results obtained with forces 
that have the same saturation properties for symmetric NM. We will refer to a "asy- 
stiff/superstiff" EOS when we are considering a potential symmetry term with a lin- 
ear/parabolic increase with nuclear density and to a "asy-soft" EOS when it shows a 
saturation and eventually a decrease above normal density 

In Figs lll2l we report, for a 124 Sn asymmetry (N — Z)/A = 0.2, the density dependence 
of the symmetry contribution to the mean-field potential (left) and of the chemical po- 
tentials (right) for neutrons (top curves) and protons (bottom curves) , for the different 
effective interactions in the isovector channel. From the behavior in the low densities 
region we expect that when the inhomogenities develop both neutrons and protons have 
the tendency to move from lower to higher density regions, in phase: the system is un- 
stable against is o scalar-like fluctuations and not isovector, see later. Since the variations 
of the two chemical potentials are different (larger for protons) we expect a lower asym- 
metry in the liquid phase. In the case of a contact between more dilute and "normal" 
density regions, we see from FigOthat in this range the neutrons have the tendency to 
move towards the dilute part producing a n-enrichment while the protons will migrate 
to the higher density regions. This mechanism is present in the "neck fragmentation", 
[ EH EJ: the neck IMF's will be always more n— rich compared to the fragments 
produced in the case of bulk fragmentation. 

Stochastic transport simulations [ [7J H2 EH E] of fragment production collisions at 
medium energies are confirming these predictions, see refs. [ EH \H\ where a comparison 
with recent data [ UHl EH is also performed. 

Isospin Distillation in Dilute Nuclear Matter 

For charge asymmetric systems we expect a qualitative new feature in the liquid-gas 
phase transition, the onset of chemical instabilities that will show up in a novel struc- 
ture of the unstable modes [ EJ- Experimentally this will be revealed through the 
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Figure 3. The mechanism of isospin distillation. The arrows show the "direction" of 
the unstable mode in various points of the spinodal region. The thick lines represent 
a composition corresponding to the initial concentration, the thin lines a pure isoscalar 
mode. Initial proton fractions y = 0.5, y = 0.4 and y = 0.1 



Isospin Fractionation or Distillation effect in multifragmentation events. We have now 
a new degree of freedom, the concentration, and in the spinodal region the fluctuations 
against which a binary system becomes unstable depend on the nature of the interaction 
between the two components of the mixture. We define density fluctuations as isoscalar- 
like in the case when proton and neutron densities fluctuate in phase and as isovector-like 
when the oscillations are out of phase. For the dilute asymmetric nuclear matter because 
of the attractive force between protons and neutrons at low density the phase transition 
is uniquely driven by isoscalar-like instabilities IT7j. 

An intuitive picture is presented in FiglHl With increasing asymmetry the direction 
of the unstable modes (arrows) in the (5p n , Sp p ) plane is more and more diverging from 
the constant concentration value (thick lines), towards a less asymmetric liquid phase. 
The angle between the two directions, i.e. the amount of isospin distillation, will be 
proportional to the repulsion of the symmetry term at sub-saturation densities. 

Symmetry effects at high baryon density: collective flows 

It is quite desirable to get information on the symmetry energy at higher density, where 
furthermore we cannot have complementary investigations from nuclear structure like in 
the case of the low density behaviour. Heavy Ion Collisions (HIC) provide a unique way 
to create asymmetric matter at high density in terrestrial laboratories. 

The isospin dependence of collective flows has been already discussed in a non-relativistic 
framework [EIIE]- The main new result shown here, in a Relativistic Mean Field (RMF) 
scheme [I19j. is the importance at higher energies of the microscopic covariant structure 
of the effective interaction in the isovector channel: effective forces with very similar 
symmetry terms can give rise to very different flows in relativistic heavy ion collisions [ 

nag. 

A full description of the isovector channel in a relativistic framework in principle should 
rely on the balance between a scalar (5 — like, attractive) and a vector (repulsive) [ 
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Figure 4. Proton-neutron differential col- 
lective flow in the 132 Sn + 132 Sn reaction 
at 1.5 AGeV b=6fm for the three differ- 
ent model for the isovector mean fields. 
Full circles and solid line: RMF — (p + 5). 
Open circles and dashed line: RMF — p. 
Stars and short dashed line : RMF — Dp. 



Figure 5. Difference between neutron and 
proton elliptic flow as a function of the 
transverse momentum in the 132 Sn+ 132 Sn 
reaction at 1.5 AGeV b=6fm in the rapid- 
ity range —0.3 < y/y pro j < 0.3. Full cir- 
cles and solid line: RMF — (p + 5). Open 
circles and dashed line: RMF — p. Stars 
and short dashed line: RMF — Dp. 



"2~1"1 123*1 IS 123 contributions. This is a quite controversial point. In relativistic HIC's, 
due to the large counterstreaming nuclear currents, one may directly exploit the different 
Lorentz nature of a scalar and a vector field [ 1320] . 

For the description of heavy ion collisions we solve the covariant transport equation 
of the Boltzmann type within the Relativistic Landau Vlasov (RLV) method [ 26J (for 
the Vlasov part) and applying a Monte-Carlo procedure for the collision term, including 
inelastic processes involving the production/absorption of nucleon resonances, [I27|. 

Typical results for the 132 Sn + 132 Sn reaction at 1.5 AGeV (semicentral collisions) are 
shown in Figs, |U [5] In FiglU we report the differential flow F pn (y) = l/N(y)J2iPx i 'Ti 
where N(y) is the total number of free nucleons at the rapidity y, p Xi is the transverse 
momentum of particle i in the reaction plane, and T\ is +1 and -1 for protons and neutrons, 
respectively. The RMF— (p+5) case (full circles and solid line) presents a stiffer behaviour 
relative to the RMF — p (open circles) model, as expected from the more repulsive 
symmetry energy E sym (p B ) at high baryon densities ["2SH2D]- We have however repeated 
the calculation using the RMF — Dp interaction, i.e. with only a p contribution but tuned 
to reproduce the same EOS of the RMF — (p + 5) case. The results, short-dashed curve of 
Figlffl are very similar to the ones of the RMF — p interaction. Therefore we can explain 
the large flow effect as mainly due to the different strengths of the vector-isovector field 
between RMF — (p + 5) and RMF — p, Dp in the relativistic dynamics. In fact if a source 
is moving the vector field is enhanced (essentially by the local 7 Lorentz factor) relative 
to the scalar one. Keeping in mind that RMF — (p + 5) has a three times larger p field it 
is clear that dynamically the vector-isovector mean field acting during the HIC is much 
greater than the one of the RMF — p, Dp cases. 

In FigJSJwe report the elliptic flow v 2 (y,pt), v 2 =< (p 2 . —p 2 y)lp 2 t > where p t = ^p 2 + p 2 
is the transverse momentum [28 . A negative value of v 2 corresponds to the emission 
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Figure 6. Central Au + Au collisions at 0.6AGeV (upper) and l.OAGeV (bottom). Time 
evolution of n/p ratio and A resonance production in high density regions (p/po > 2.0) 
(first two columns) and of the total ir~ /ir + ratio (right). Solid lines: RMF — p. Dashed 
lines: RMF - (p + 5). 



of matter perpendicular to the reaction plane, sqeeze — out flow. The ^-dependence 
of t>2 is very sensitive to the high density behavior of the EOS since highly energetic 
particles (jp t > 0.5) originate from the initial compressed and out-of-equilibrium phase of 
the collision. We focus on the proton-neutron difference of the elliptic flow. From Fig|S]we 
see that in the (p + 5) dynamics the high-p t neutrons show a much larger squeeze — out. 
This is fully consistent with an early emission (more spectator shadowing) due to the 
larger repulsive p-field. The v 2 observable, which is a good chronometer of the reaction 
dynamics, appears to be particularly sensitive to the Lorentz structure of the effective 
interaction. 
7t~/tt + ratios 

Using the same relativistic transport code we have evaluated the tc~ vs. tt + production 
for central Au + Au collisions at different energies, see FigEl As expected the larger 
repulsion seen by neutrons at high densities in the RMF — (p + 5) will show up in a 
reduced n/p ratio, smaller A ' - density and finally a reduced 7r~ production. However 
now reabsorption effects are important and actually the effect appears to be decreasing 
at higher energies. 

Isospin and Deconfinement at High Baryon Density 

It is relatively easy to show that at high baryon density and low temperature we can 
expect a transition from hadronic matter to deconfined quark matter. The procedure is 
straightforward: 

(i) Start from two "reasonable" model Equations of State (EOS), one for the hadronic 
phase, which can reproduce saturation properties, one for the quark phase, which can 
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Figure 7. Variation of the transition density with proton fraction followed in the interac- 
tion zone during a semicentral 132 Sn + 132 Sn collision at lAGeV (circles) and SOOAMeV 
(crosses) 



reproduce the hadron spectrum. 

(ii) Construct the phase separation boundary surface from the Gibbs phase rule. 

For symmetric matter the baryon density p tr corresponding to the transition to the 
coexistence region is relatively high, as expected, ranging from 4 to 8 times the saturation 
value po, depending on the stiffness of the hadronic EOS at high densities. The new 
feature we would like to focus on in this report is the isospin dependence of such boundary 
location. We can foresee an interesting asymmetry effect, in the appealing direction of a 
decrease of p tr , since the hadronic EOS becomes more repulsive ref . [ [2H1 EH] • 

The proton fraction Z/A dependence of the p tr is reported in Fig [7] with the bag constant 
value B l l A = 150 MeV and a s = for the quark phase and various choices for the hadronic 
EOS: Dotted line GM3 parametrization [ 31J; Dashed line RMF — p parametrization [ 
123] : Solid line RMF — (p + 5) parametrization GM3 and RMF — p have the same 

source of the interaction symmetry term (only the p-meson). The effects of the asymmetry 
appears now quite dramatic: we see a ptrans as low as 2po for proton fractions between 
0.3 and 0.4, conditions that with some confidence we could "locally" reach in a heavy ion 
collision at intermediate energy using exotic very asymmetric beams. 

Using our Relativistic Transport Code, with the RMF — (p + 5) effective interaction, 
we have performed some simulations of the 132 Sn + 132 Sn (average Z/A = 0.38) collision 
at various energies, for a semicentral impact parameter, b = 6/m, just to optimize 
the neutron skin effect. In Fig|7| the paths in the (p,Z/A) plane followed in the cm. 
region during the collision are reported, at two energies 300 AMeV (crosses) and 1 AGeV 
(circles). We see that already at 300 AMeV we are reaching the border of the mixed 
phase, and we are well inside it at 1 AGeV . 

In conclusion we support the possibility of observing precursor signals of the phase 
transition to a deconfined matter in violent collision (central and semicentral) of exotic 
(radioactive) heavy ions in the energy range of few hundred MeV per nucleon. A possible 
signature could be revealed through an earlier " softening" of the hadronic EOS for larger 
asymmetries. 
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Outlook: The Eleven Observables 

As a conclusion of our report we like to suggest a selection of Eleven Observables, 
from low to relativistic energies, that we expect particularly sensitive to the microscopic 
structure of the in medium interaction in the isovector channel, i.e. to the symmetry 
energy and its "fine structure": 

1. Competition of Reaction Mechanisms. Interplay of low-energy dissipative mechanisms, 
e.g. [ 7, fusion (incomplete) vs. deep-inelastic vs. neck fragmentation: a stiff symmetry 
term leads to a more repulsive dynamics. 

2. N/Z of fast nucleon emission. Symmetry repulsion of the neutron/proton mean field 
in various density regions. 

3. Neutron/ Proton correlation functions. Time-space structure of the fast particle emis- 
sion and its relation to the baryon density of the source, see the recent [ I32j . 

4- Fragment Multiplicities. A more efficient use of protons in forming primary fragments 
is expected in the asy-stiff case. 

5. Isospin Distillation (Fractionation). Isospin content of the Intermediate Mass Frag- 
ments in central collisions. Test of the symmetry term in dilute matter. 

6. Isospin content of Neck- Fragments. Test of the symmetry term around po- 

7. Fast Fission Multiplicity. The rate of "aligned" fission events of the Projectile- 
Like/Target-Like Fragments reflects the symmetry repulsion in semicentral collisions. 

8. Isospin Diffusion. Measure of charge equilibration in the "spectator" region in semi- 
central collisions, test of symmetry repulsion. 

9. Neutron-Proton Collective Flows. Together with light isobar flows. Check of symme- 
try transport effects. Test of the momentum dependence (relativistic structure) of the 
interaction in the isovector channel. Measurements also for different p t selections. 

10. 7i~ / 7r + Yields. Since ir~ are mostly produced in nn collisions we can expect a reduction 
for highly repulsive symmetry terms at high baryon density, see [ £21 EH] • 

11. Deconfinement Precursors. Signals of a mixed phase formation (quark-bubbles) in 
high baryon density regions reached with asymmetric HIC at intermediate energies. 

From the points 3, 4, 5, 8, 9 in our simulations we presently get some indications for asy- 
stiff behaviors, i.e. increasing repulsive density dependence of the symmetry term, but 
not more fundamental details. Moreover all the available data are obtained with stable 
beams, i.e. within low asymmetries. 
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